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Abstract. Very- high energy photons emitted by distant cosmic sources are absorbed 
on the extragalactic background light (EBL) during their propagation. This effect can 
be characterized in terms of a photon transfer function at Earth. The presence of 
extragalactic magnetic fields could also induce conversions between very high-energy 
photons and hypothetical axion-like particles (ALPs). The turbulent structure of the 
extragalactic magnetic fields would produce a stochastic behaviour in these conversions, 
leading to a statistical distribution of the photon transfer functions for the different 
realizations of the random magnetic fields. To characterize this effect, we derive new 
equations to calculate the mean and the variance of this distribution. We find that, 
in presence of ALP conversions, the photon transfer functions on different lines of 
sight could have relevant deviations with respect to the mean value, producing both 
an enhancement or a suppression in the observable photon flux with respect to the 
expectations with only absorption. As a consequence, the most striking signature 
of the mixing with ALPs would be a reconstructed EBL density from TeV photon 
observations which appears to vary over different directions of the sky: consistent with 
standard expectations in some regions, but inconsistent in others. 
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1. Introduction 

Axion-like particles (ALPs) with a two-photon vertex are predicted in many extensions 
of the Standard Model [1-3]. The CZ77 coupling allows for ALP-photon conversions in 
electric or magnetic field. This effect is exploited by the ADMX experiment to search 
for axion dark matter [4], by CAST to search for solar axions [5-7], and by regeneration 
laser experiments [8-12]. 

ALPs also play an intriguing role in astrophysics. Indeed, photons emitted by 
distant sources and propagating through cosmic magnetic fields can oscillate into ALPs. 
The consequences of this effect have been studied in different situations [13-23]. In 
particular, in the last recent years photon-ALP conversions have been proposed as a 
mechanism to avoid the opacity of the extragalactic sky to high-energy radiation due to 
pair production on the Extragalactic Background Light (EBL). At this regard, recent 
observations of cosmologically distant gamma-ray sources by ground-based gamma-ray 
telescopes have revealed a surprising degree of transparency of the universe to very high- 
energy (VHE) photons (E > 100 GeV) [24, 25]. Surprisingly, data seem to require a lower 
density of the EBL than expected and/or considerably harder injection spectra than 
initially thought [26, 27]. Oscillations between very high-energy photons and ALPs could 
represent an intriguing possibility to explain this puzzle through a sort of "cosmic light- 
shining through wall" effect. Infact, if VHE photons are converted into ALPS and then 
regenerated, they should not suffer absorption effects while they propagate as ALPs. In 
this sense, two complementary mechanisms have been proposed: a) VHE photon-ALP 
conversions in the magnetic fields around gamma-ray sources [20, 28] and then further 
back-conversions in the magnetic field of the Milky Way [29] (see also [30]); b) oscillations 
of VHE photons into ALPs in the random extragalactic magnetic fields [31,32]. In 
principle, both the mechanisms can be combined together, as shown in [33]. Currently, 
the inference of EBL from VHE photons emitted by sources at high redshift [34] is still 
object of debate, and it is not clear how robust are the conclusions on the absorption 
effects obtained from the recent gamma data [35]. In this sense, it is also possible that 
the observed transparency of the universe to VHE photons could be explained without 
the need of introducing nonstandard mechanisms (see, e.g., [36,37]). Nevertheless, the 
seminal works mentioned before have pointed out the nice connection between VHE 
gamma-astronomy and ALP searches. Therefore, it seems worthwhile to further explore 
the consequences of this exciting possibility. 

In this context, the treatment of the oscillations of VHE photons into ALPs in 
the extragalactic medium in presence of the absorption on the EBL presents a certain 
degree of complexity. Extragalactic magnetic fields are supposed to have a turbulent 
structure which can dramatically affect the development of photon-ALP conversions. 
Brute force numerical simulations get the solution of the mixing equations along a given 
photon line of sight by iterating the equations in each domain in which the magnetic 
field is assumed constant [14,31]. Since one cannot know the given configuration of 
magnetic domains crossed by VHE photons during their propagation, in the previous 
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literature VHE photon-ALP conversions were usually characterized in terms of the 
mean conversion probability, obtained averaging the resulting conversion probabilities 
over an ensemble of magnetic field configurations along the photon line of sight. Such 
a procedure slows down the solution of our problem, since to obtain stable results 
typically the average has to be performed over more than 10 3 realizations of the magnetic 
fields [31]. Moreover, the use of the mean probability as representative value does not 
appear a priori completely justified, since the variance of the probability distribution 
could produce relevant deviations from the mean value for conversions occurring in 
different realizations of the random magnetic fields. 

In order to overcome these previous limitations, in this paper we perform a new 
study of the mixing equations of VHE photons in the turbulent magnetic fields and we 
provide a user-friendly calculation of the mean and of the variance for the distribution 
of the photon transfer functions. The plan of our work is as follows. In Section 2 we 
discuss about the absorption of VHE photons on the extragalactic background light. In 
Section 3 we characterize VHE photon-ALP mixing in presence of absorption and we 
present our calculation of the mean photon transfer function averaged over the ensemble 
of all the possible realizations of random magnetic field configurations. We also show 
how to calculate the variance of the statistical distribution of the transfer functions. In 
Section 4 we present our results for the photon transfer function of VHE gamma-rays 
emitted from distant sources with and without the mixing with ALPs. Contrarily to 
previous predictions, the presence of a broad variance in the statistical distribution of the 
photon transfer functions could produce both an enhancement or a suppression of the 
observed VHE photon flux with respect to the case with only absorption. The resulting 
photon flux would depend on the particular random magnetic field configuration along 
the photon line of sight. As a consequence, photon-ALP mixing can not provide an 
universal mechanism to obtain the transparency of the universe to VHE radiation, but 
instead they would produce a strong direction-dependent behaviour in the flux of VHE 
photons from distant sources. In Section 5 we discuss about possible developments of 
our study and we conclude. There follow two Appendices, in which we present some 
details for the derivation of the mean photon transfer function (Appendix A) and for 
the variance of the distribution (Appendix B). 

2. Absorption of very high-energy photons on extragalactic background 
light 

The flux of very high-energy (VHE) gamma rays (E > 100 GeV) from distant sources 
is attenuated in an energy dependent way by the interaction with background photons 
in the universe. The main source of absorption for VHE photons is due to the pair 
production process ^ VHE ^ bk s —> e + e~. In the energy range 100 GeV <E< 10 TeV, 
the absorption is dominated by the interaction with optical/infrared photons of the so 
called Extragalactic Background Light (EBL), sometimes also referred as Metagalactic 
Radiation Field (MRF). The absorption rate for such a process in function of the incident 
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photon energy E is given by [38] 
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where the limits of integration in both integrals are determined by the kinematical 
threshold of the process and 
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with o"o = 1.25 x 10 -25 cm 2 , is the cross section for the pair production process [39], 
in function of the electron velocity in the center of mass of the interaction (3 = 
[1 — 2m 2 /£'e(l — ^)} 1 ^ 2 , being e the background photon energy, and £ the cosine of 
the angle between the incident and the background photon. For practical purposes, it 
can be useful to notice that the inner integral in d£ in Eq. (1) can be evaluated by 
performing the change of variable £ = £(/?) and has actually an analytic closed form: 
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where (3 m = (1 — m 2 e / Ee) 1 ^ 2 is the maximum electron velocity in the center of mass of 
the interaction and the function Li 2 is the polylogarithm of order two. 

For simple estimations, the background photon spectrum in Eq. (1) can be 
approximated, at redshift z = 0, with a power-law [40] 

7 -10-^(47) eV^cm" 3 , (4) 




de VeV, 

with 0.61 < k < 1.52, depending on the model used (see Fig. 1). With this 
approximation one has 

TJE) / E \ L55 

7V i ~ 1.1 x 10~ 3 fc — — . 5 
Mpc" 1 VTeVy v ; 

In the literature are present different realistic models for the photon background. In 
particular, in the following we will refer to the 2008 Minimal Kneiske Model [41], which 
provides a strict lower-limit flux for the extragalactic background light from ultraviolet 
to the far-infrared photon energies. The model parameters are chosen to fit the lower 
limit data from galaxy number count observations, assuming that the shape and the 
normalization of the background radiation does not change except for the red-shifting 
between the source and the observer. This model has the advantage that it is not inferred 
by an inversion of the VHE photon observed spectra, as e.g. in [25]. We note that this 
latter procedure would be not reliable in the presence of VHE photon-ALP conversions, 
since in this case the effects of mixing and absorption on EBL would be entangled, 
preventing the possibility to extract information on the EBL from VHE photon measured 
spectra. Our reference model gives us the maximal possible transparency compatible 
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Figure 1. Spectrum of the EBL relevant for the absorption of VHE photons. The 
dashed power-law lines correspond to the simplified model of Eq. (4) for the limiting 
cases of k = 0.61 (lower line) and k = 1.52 (upper line) respectively. The solid line 
corresponds to the 2008 Minimal Kneiske Model at rcdshift z = 0. (see the text for 
details) 



with the standard expectations, so that an evidence of a greater transparency would 
have to be attributed to nonstandard effects in the photon propagation. 

The spectral energy distribution at a given redshift z can be inferred by the 
tabulated power spectrum [41] P(X, z) = z) (where /(A, z) is the flux at redshift 

z of energy for unit of solid angle between A and A + d\ where A is the comoving 
wavelength) by 

d ^ = 2P{\{l + z)){l + zf. (6) 

The corresponding EBL energy spectrum is represented in Fig. 1. From this Figure 
we see that the realistic Kneiske model has an approximate power-law trend at photon 
background energies e> few xlO -1 eV, which are relevant for the absorption of TeV 
photons. This will be useful to explain the high-energy behaviour of the transfer 
function. 
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The absorption function T 7 as function of the photon energy on Earth E and of 
the redshift z is given by Eq. (1) simply changing E — > E(l + z) in the integrand and 
in the limit of integration. 

Finally, for a given source at distance L, the photon spectrum observed on Earth 
(apart from the geometrical dilution) is given by 

Iobs(E) = T 7 (E, L) ■ / S ourcc(-Eo) = exp (-T 7 ) J SO urce(-Eo) , (7) 

where the initial source spectrum /source (-^0)) with initial photon energy E = E(l+z), is 

modified due to the effect of the VHE photon absorption, through the transfer function 

T 7 , expressed in terms of the optical depth 

rL r r z o dz 
r 7 =/ dxT 1 (E,x) = — T 7 (E,z), (8) 

Jo H * Jo (i + z)Jn A + n m n + z) 3 



dependent on the evolution of the Universe through the Hubble constant Hq = 
73 km Mpc -1 s _1 , the matter density Q m = 0.24 [42], and the dark energy density 
Q\ = 1 — Q m (assuming a flat cosmology). We observe that in absence of photon- ALP 
conversions the transfer function T 7 drops exponentially at high energies. We will see 
that the effect of VHE photon- ALP conversion is to soften this behaviour. 



3. Very high-energy photons mixing with axion-like particles 

3.1. Equations of motion 

The effect of the absorption of VHE photons on the EBL, described in the previous 
Section, can be strongly modified if photons do mix with axion-like particles (ALPs). 
Pseudoscalar ALPs couple with photons through the following effective Lagrangian [43] 
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(9) 



where a is the ALP field with mass m a , F^ the electromagnetic field-strength tensor, 
F^y = \e^ upa F p<7 its dual, and g ai the ALP-photon coupling. As a consequence of this 
coupling, ALPs and photons do oscillate into each other in an external magnetic field. 
For a scalar particle, the coupling is proportional to F^F^a. For definiteness, we limit 
our discussion to the pseudoscalar case, but similar consequences apply also to scalars. 

Let us suppose that a photon with energy E moves in the x% direction. The 
transverse component of the external magnetic field is = B — B^e^. The evolution 
equations of the photon-ALP system in presence of mixing and absorption are [14, 43] 
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where = cos0 = By • Bx/Bt = J 1 — s^. The entries A^- (i, j = 1,2) that mix the 
photon polarization states are energy-dependent terms determined by the properties 
of the medium and the QED vacuum polarization effect. In particular, neglecting the 
Faraday rotation effects which are not relevant for the high energies of our interest, they 
read 



A n = A||4 + A ± 4, 
A 22 = A|,4 + A ± 4, 

A 12 = A 21 = (Ay - A ± ) w , (11) 
with Ay = A p i + f A QED , Aj_ = A pl + 2A QED and 

A a7 = -g ai B T ~ 1.52 x 10" 2 ( T ) ( -^—) Mpc" 1 , 

7 2 ya7 VlO^GeV" 1 / V10" 9 G/ P ' 

™% d ( m a \ 2 ( E \ ~ l , 

^ni n ( E \ 1 f n e 



2 



A -^(^J— °- 9 (4)(t^gJ^-'' < 12 > 

where -Bt is expressed in Lorentz-Heaviside units and = A7ran e /m e is the plasma 
frequency of the medium, being n e the electron density. For the numerical estimations 
above, that we will use in the following as benchmark values, we have referred to 
the following physical input: The strength of widespread, all-pervading 5-fields in 
the extragalactic medium must be B <2.8 x 10~ 7 (//Mpc) _1/2 G, coherent on a scale 
I ~ 1 Mpc [44], as obtained scaling the original bound from the Faraday effect of distant 
radio sources [45, 46] to the now much better known baryon density measured by the 
Wilkinson Microwave Anisotropy Probe (WMAP) [47]. The mean diffuse intergalactic 
plasma density is bounded by n e <^ 2.7 x 10~ 7 cm -3 , corresponding to the recent WMAP 
measurement of the baryon density [47]. Recent results from the CAST experiment give 
a direct experimental bound on the ALP-photon coupling of g^ £8.8 x 1CT 11 GeV" 1 
for m a <0.02 eV [7], slightly better than the long-standing globular-cluster limit [48]. 
For ultra-light axions a stringent limit from the absence of 7-rays from SN 1987A gives 
g ai £ 1 x 1CT 11 GeV" 1 [49] or even g aj £ 3 x 1CT 12 GeV" 1 [50]. Previous bounds on ALPs 
can be relaxed if they have a chameleontic nature [51]. In this case, the best constraint 
comes from the structure of starlight polarization: g ai <> 10" 9 GeV" 1 [52]. 

The absorption term T 7 in Eq. (10), due to the VHE photons scattering with the 
low-energy photons in the background, produces a damping of the oscillations in analogy 
with the case of the mixing of high-energy neutrinos in an absorbing matter [53] . In the 
presence of this term, the Hamiltonian 7i is no longer hermitian. 
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3.2. Mean photon transfer function 



Assuming an homogeneous magnetic field in a domain of size /, in absence of absorption 
the probability that a photon will convert into an ALP reads [43] 
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2A a7 , the photon-ALP mixing is maximal 
[9 ~ tt/4) and the conversion probability becomes energy- independent. This is the so- 
called strong- mixing regime. In this case, if A aj l <C 1, the conversion probability on a 
single domain becomes extremely simple, namely P ai = (A a7 /) 2 . In the following, we 
will work in this regime. 

In this situation, we explicitly drop Ay^ and A a from the equations of motion 
[Eq. (10)]. Thus, the propagation hamiltonian 7i can be written as 7i = A — iD where 
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and V = ^diag(l, 1, 0) is the damping term associated to the absorption. 

Very high-energy gamma-rays propagate in the extragalactic magnetic fields during 
their route to the Earth. These £>-fields presumably have a turbulent structure. 
Therefore, for the case under study we need to describe photon-ALP conversions in 
random magnetic field configurations. Let us now consider the propagation of photons 
in many domains of equal size I (~ 1 Mpc in our case) in which the magnetic field has 
(constant) random values and directions. Along a given line of sight, the angles are 
randomly distributed in [0, 2ir). In the following, we will work in the formalism of the 
density matrix 

Mi \ 
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For the k-ih domain the density matrix is given by 

Pk = e~ iUkl ■ Pk-i ■ e m l l , (19) 

where e~ tHkl is the propagation operator for the fc-th domain. During their path with 
a total length L, photons cross k = l,...n domains (n = L/l) representing a given 
random realization of and 0*.. Since we cannot know this particular configuration, 
we perform an ensemble average over all the possible realizations on the 1, . . . n domains. 
Defining this ensemble average as p n = (p n )i... n , we have 

Pn = (e~ m " 1 ■ Pn-i ■ e m " l h.. n = (e- m " 1 ■ p n ., ■ e»&) n . (20) 

For the chosen values of the input parameters as in Eq. (12), we can perform a 
perturbative expansion up to the second order of the evolution operator in each domain, 
i.e. 

e -m n i _ 1 _ m j _ }_ H 2j2 _ ( 21 ) 

Performing then the ensemble average, as shown in Appendix A, using that 
Pn — Pn-i — l9 X3 p(xs), and summing over the two indistinguishable photon polarization 
states, we finally arrive at a system of two coupled differential equations 



d ( T \ P 
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(22) 



where T 7 = pn + p 2 2 and T a = p aa are the mean transfer functions for the photon and 
for the ALP respectively; P ai = A^ 7 / 2 is the average photon- ALP conversion probability 
in each domain (in absence of absorption and in the limit of strong mixing) and 
finally a = r 7 //P a7 is the ratio between the absorption probability and the conversion 
probability. 

In realistic astrophysical situations both P ai and a are functions of the distance, 
due to the redshift dependence of the extragalactic magnetic field and of the EBL. 
However, taking these parameters as constant, in the hypothesis of only photons in the 
initial state (T 7 (0) = 1, T a (0) = 0) Eq. (22) has a simple analytical solution 



T 7 (j/) = e" 
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where 
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K = V V 2 — OL , 

y = . (24) 
In particular, Eq. (23) gives the two limiting expressions 

f 2 , 1 -3y/2 o 

a »;. < 25 ' 

In the absence of absorption, we recover the mean transfer function already found in 
[54]. Conversely, in the case of strong absorption (a ^> 1) we have T 7 oc (r 7 )~ 2 . Using 
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Figure 2. Photon transfer function T 7 (y) for a = (no absorption) in the upper 
panels and for a = 1 in the lower panels. The continuous curve corresponds to the 
mean value obtained by Eq. (23), while the dashed lines correspond to the dispersion 
around the mean T 7 ±(5T 7 , calculated in Appendix B. In the left panels, the dot-dashed 
lines correspond to a given realization of the random magnetic fields along the photon 
line of sight, obtained by the numerical integration of Eq. (10). In the right panels 
we show a scatter plot for T 7 corresponding to N r = 20 realizations of the random 
magnetic fields, (see the text for details) 



the approximate expression for T 7 given in Eq. (4) we observe that the transfer function 
would drop as a power of the energy (rather than exponentially as expected without 
ALP mixing). Moreover, also the attenuation of the transfer function with the distance 
is less than in the case of absence of conversions. In fact, the argument of the exponential 
is suppressed by a factor 1/at with respect to the no-conversion case [see Eq. (7)]. 

In Appendix B we report also the calculation of the root mean square <5T 7 for the 
distribution of the transfer function in different random realizations of the magnetic 
field. This result is useful to estimate the uncertainty associated with the averaging 
procedure. 

For illustrative purposes, in Fig. 2 we compare the transfer matrix T 7 (?/) of Eq. (23), 
with the numerical solution of Eq. (10), for a value of the parameter A a7 as in Eq. (12). 
For simplicity, we choose constant values of the absorption factor a, namely a = 
(no absorption) in the upper panel and a = 1 in the lower panel. The dashed lines 
represent T 7 ± <5T 7 , where <5T 7 is the dispersion calculated as in Appendix B. In the 
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left panels, the dot-dashed lines correspond to T 7 for a given random realization of 
the magnetic field along the photon line of sight, obtained by the numerical solution 
of Eq. (10). In absence of absorption, we realize that along a given line of sight, the 
photon transfer function can present strong deviations with respect to the average value 
obtained with our analytical calculation. In this case, the dispersion with respect to 
the average tends to 5T 7 = l/3y/E for y — > oo, as shown in Appendix B. We stress that 
the presence of a dispersion around the average T 7 has not been properly appreciated 
in the previous literature. In particular, previous studies on the mixing of photons with 
ALPs emitted from point-like sources into random magnetic fields have presented 2/3 
as limiting value for the photon transfer function (see e.g. [13] for the case of photons 
emitted by supernovae la). In the case of strong mixing, where the oscillations are 
achromatic and photons of different energies are not dephased, this result is correct 
only on average, while along a given line of sight one can expect 0(1) deviations from 
this value. This peculiar effect has been also recently recognized in the context of 
photons emitted by active galactic nuclei [23]. Indeed, the scatter in their luminosities 
has been interpreted in [23] as a possible hint of the existence of a very light ALP. As a 
further consequence of this effect, it would be worthwhile to investigate if the presence 
of these large variations in the photon transfer function could put additional constraints 
to the mechanism of photon- ALP conversions introduced in [13] to explain the observed 
dimming of supernovae la. In fact, photon conversions into ALP could produce large 
dispersions in the supernova lightcurves. In the presence of absorption (lower panels) 
the transfer function (and its dispersion) is suppressed as a power- law for y ^> 1, due to 
the damping of the oscillations. In the right panels, we superimpose on our prescription 
for T 7 a scatter plot corresponding to N r = 20 different realizations of the random 
magnetic fields and for 500 steps in the variable y. Again, we realize that behaviour 
of the T 7 distribution is damped by the effect of the absorption (a = 1). However, as 
we will see in the next Section, the presence of a dispersion around the mean value will 
have interesting consequences for VHE photons. Finally, we mention that performing 
the averaging procedure over the transfer functions corresponding the different random 
realizations, we recover our analytical results (not shown). 

4. Transfer function for very high-energy photons 

We will apply the results of the previous Section to the case of VHE photons. For the 
characterization of the EBL, we refer to the 2008 Minimal Kneiske Model [41], discussed 
in Section 2. Concerning the extragalactic magnetic field, we consider it frozen into the 
medium. With this assumption, the scaling law of the magnetic field with the redshift 
z is given by B(z) = Bq(1 + z) 2 [46] while the size of the magnetic domains scales as 
I = /o/(l + z). With this choice, P a - t {z) scales as (1 + z) 2 . We assume B = 1 nG and 
l — 1 Mpc. Equation (22) can be easily written in terms of the redshift, by means of 
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In Fig. 3 we show the photon transfer function T 7 in presence of absorption on the 
EBL for our reference model, in function of the redshift z for different values of the 
observed photon energy E with (continuous curve) and without ALP mixing (dotted 
curve). The dashed curves represent the spread T 7 ± <5T 7 around the mean value. 
One realizes that at redshift z > 0.2 the presence of ALPs could produce dramatic 
modifications in the shape of the photon transfer function, the stronger the effect the 
higher the photon energy. In this sense, a hint of the mixing of very high-energy photons 
with ALPs would be the detection of very distant sources otherwise obscure due to the 
absorption. However, from the spread of the photon transfer function <5T 7 , we also realize 
that at high redshift T 7 can present relevant deviations with respect to the mean value. 
In this sense, the effect of mixing with ALPs for VHE photons emitted by distant gamma 
sources would be strongly dependent on the particular realization of the extragalactic 
magnetic fields crossed by them during their propagation. Therefore, it is not guaranteed 
that the mixing with ALPs could provide an universal mechanism to achieve the 
transparency of distant gamma-sources. Conversely, for particular configurations of the 
extragalactic magnetic fields crossed by VHE photons, one could also find a suppression 
of the transfer function stronger than in the standard case. In general, observing VHE 
photons from far sources one would find strong differences in the measured spectra along 
different lines of sight. As a further consequence, the presence of ALPs would prevent the 
possibility to use distant gamma sources at very high energy as cosmological candles [55]. 

In Fig. 4 we show the photon transfer function as function of the energy for four 
different values of the redshift of the emitting source (as in Fig. 1 of [41]). We see that in 
absence of ALP conversions the photon transfer function would be strongly suppressed 
at energies above E > 100 GeV, the stronger the suppression the larger the redshift. On 
the other hand, in presence of conversions, T 7 has an approximate power-law behaviour 
at high energies, as pointed out in the previous Section. We also realize that the spread 
in the possible values of T 7 would make difficult to infer strong conclusions about ALP 
mixing observing only few sources. We also note that in the case of z = 0.20 the inclusion 
of the ALPs does not produce any significant change in the photon transfer function. 
This suggests that it would be difficult to interpret in terms of ALP conversions the 
presumed transparency to gamma radiations for the sources at z = 0.165 and z = 0.186 
discussed in [24,25]. Conversely, ALP conversions could play a significant role for the 
source 3C279 at redshift z = 0.54 [34]. 

Finally, in Figure 5 we show two iso-contours of the photon transfer function T 7 in 
function of the redshift z and of the observed photon energy E with and without ALP 
conversions. We see that the attenuation of the photon transfer function in the presence 
of ALP conversions could present a large variation at high redshift with respect to the 
standard expectations with only absorption. 
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Figure 3. VHE photons transfer function T 7 in function of the redshift z, for different 
values of the observed photon energy E. The standard case, including only photon 
absorption on the EBL, is represented by the dot-dashed curve. In the presence 
of mixing with ALPs the mean photon transfer function, obtained by Eq. (22), is 
represented by the solid curve and its typical spread T 7 ± <5T 7 corresponds to the two 
dashed lines. 



5. Discussion and conclusions 

Very high-energy gamma-ray observations would open the possibility to probe the 
existence of axion-like particles (ALPs) predicted in many theories beyond the Standard 
Model. Recent gamma observations of cosmologically distant gamma-ray sources 
have revealed a surprising degree of transparency of the universe to very high-energy 
photons. The oscillations between high-energy photons and ALPs in the random 
extragalactic magnetic fields have been proposed as an intriguing possibility to explain 
these observations [31,32]. Apart from the original proposal, the consequences of 
this mechanism are testable with the measurements of the new generation of Imaging 
Atmospheric Cherenkov Telescopes, like MAGIC [56], HESS [57], VERITAS [58] or 
CANGAROO-III [59], covering energies in the range 0.1-20 TeV, and hopefully with 
the future Cherenkov Telescope Array, reaching energies of 100 TeV [60]. 

In order to perform a systematic study of these effects, without recurring to brute- 
force time-consuming numerical simulations, in this paper we have presented a simple 
calculations of the photon mean transfer function and of its variance in presence of 
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log 10 (E/ GeV) log 10 (E/ GeV) 

Figure 4. VHE photons transfer function T 7 in function of the observed photon 
energy E, for different values of the redshift z. The labels of the different curves are 
the same as in Fig. 3. 

absorption on the EBL and mixing with ALPs. We have found that our prescription is 
enough accurate for the case under study and we have shown some numerical examples 
of VHE photon transfer functions. It results that VHE photon-ALP mixing would 
produce peculiar deformations in the energy spectra of very high-energy gamma-rays 
emitted at high redshift. We have also found that the photon transfer function in 
presence of VHE photon-ALP conversions presents a relevant dispersion around the 
mean value due to the randomness of the extragalactic magnetic fields crossed by the 
photons. Due to this fact, the measured flux for VHE photons traveling along different 
lines of sight can be strongly suppressed or enhanced with respect to the case with 
only absorption, depending on the particular configuration of random magnetic fields 
encountered. This would suggest that photon-ALP conversions could not be an universal 
mechanism to produce the transparency of universe to VHE photons. Conversely, the 
most striking signature of the mixing with ALPs would be a reconstructed EBL density 
from TeV photon observations which appears to vary over different directions of the sky: 
consistent with standard expectations in some regions, but inconsistent in others. To 
test this effect we would need to collect data from sources along different directions in the 
sky in order to perform a study of the photon energy distributions, from which we could 
hope to infer possible hints of ALPs. A further signature of these stochastic conversions 
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Figure 5. Iso-contours of the photon transfer function T 7 in function of the redshift z 
and of the observed photon energy with only absorption (light curves) and in presence 
of ALP conversions. In this latter case, we represent the mean value of T 7 (black solid 
curves) and its dispersion T 7 ± 5T 7 (black dashed curves). 



would be the detection of peculiar direction-dependent dimming effects in the diffuse 
photon radiation observable in GeV range, testable with the FERMI (previously called 
GLAST) experiment [61]. 

As further developments, we plan to use our calculation to perform a systematic 
study of ALP signatures in very high-energy gamma-rays, analyzing in details the 
spectral deformations expected for observed sources at different redshifts and for 
different models of the extragalactic background light. Thanks to our calculation, this 
task now appears more doable than before. After that, it will remain to see if elusive 
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ALPs will show up from the sky. 
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Appendix A. Calculation of the mean photon transfer function 

Here, we present the derivation of the photon transfer function T 7 introduced in Sec. 3.2. 
We have defined p n = (p n ) ensemble average of the density matrix p n in the 

n-th domain over all the possible realizations of the random magnetic fields in all the 
domains from 1 to n. Since 7i n depends only from the configuration of the n-th domain, 
we have that 



Pn = (e- lHnl ■ Pn-i ■ e m " L ) 



l...n 



(e 



-iTi.nl 



Pn-i ■ e*™). 



(a.i; 



For the chosen values of the input parameters in Eq. (10), we can perform a perturbative 
expansion of the evolution operator up to the second order 7i n = A,„ — %D n in each 
domain, i.e. 
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Since, from Eq. (17) it results that (A n )0 n = 0, Eq. (A.I) can be written as 

Pn = Pn-l - I (DnPn-1 + P~n-<Pn) 

+ / 2 (A„p n _iA n )„ - - ((A 2 ) n p n _i + Pn-i(A^), 

In the previous equation we have neglected the second order terms in T> n l since they give 
a contribution of the order of (r 7 /) 2 ^ 10~ 6 [see Eq. (5)] which is at least two orders of 
magnitude smaller than (A ao i) 2 ~ 10~ 4 , when we use the benchmark values in Eq. (12). 
For (A 2 ) n we have 
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(A.4) 
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where we have defined 

1 
4 



A 2 



1 ^ 2 



(A.5) 



where _B 2 ff 



(•S|>)b„ = 2(|B| 2 )/3 due to the projection effect. In the same way we have 



(A n p n _iA„) n = -A 2 7 



Paa p a i 
Paa Pa2 
Pla P2a Pll + P22 

After n domains, the total distance travelled by the photon is x 3 = nl. Defining 
p n = p(x 3 ), one has p n — p n _i ~ W X3 p(x 3 ). After a straightforward calculation, we 
obtain the evolution equation for the averaged density matrix p 
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where 



H = a + 
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and we have defined as P. 



u-7 



A 2 7 / 2 the average photon-ALP oscillation probability 
in each domain (in absence of absorption and in the limit of strong mixing) and 
a = V y l/P ai the ratio between the absorption probability and the oscillation probability. 
Since we are interested in determining the total final photon and ALP flux, we 



define the mean transfer functions T 7 
finally obtain 
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Defining the variable dy = P aj dx 3 /l, we can also obtain a second order equation for the 



function 71 
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where the prime denotes the derivative with respect to y. For constant a we easily 
obtain the solution in Eq. (23). 



Appendix B. Calculation of higher momenta 

In the previous Appendix we have have calculated the photon transfer function T 7 
averaged over all the possible magnetic domain configurations. However, photons 
coming from a single source cross just one (unknown) particular realization of the 
magnetic field domains. It is thus interesting to evaluate the uncertainty introduced 
by the procedure of averaging. To do this, we calculate the second order momenta of 
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the probability distributions, by means of a procedure similar to the one introduced in 
Appendix A. We define the "square" of the density matrix as 

p( 2 ) = p p _> p^l = PijPkl , (B. 1) 

and p~W = (pffl)i.., n - We rewrite Eq. (A. 7) in the following way 

d _ 

where the tensor can be written as 
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otherwise . 



Performing the average of p^ up to the first order in T 7 Z and to the second order in 
A ai l as in Eq. (A. 3), we arrive at the following equation for the evolution of p^ 
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where for simplicity we have dropped the subscript n from the averages in the last 
two lines. The last term arises from the linear term in I in Eq. (20) which is 
absent in Eq. (A. 3) since it averages out ((A) = 0). From Eq. (17) we have that 
(AyA w ) = l^iaki with 

!1 if ijkl = llaa, laal, alia, aall 
22aa, 2aa2, a22a, aa22 
otherwise . 

After a long but straightforward derivation, one can extract a subset of 6 independent 
equations out from the set of the 81 of (B.3): 
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where, as usual, dy = P ai dx 3 /l and, 

R-y = Pirn + P2222 + P1122 + P2211 = ((Pn + P22) 2 ) 
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We can thus give a physical interpretation for some of these quantities: i? 7 and R a 
are the square average of the photon and ALP transfer function respectively, R p is the 
average square degree of polarization of the photon, r/ a7 is the photon-ALP correlation. 
Since the third and the fourth of Eqs. (B.4) are similar, starting from a completely 
unpolarized photon state we find also that R p (y) = (-y(y), so that we can reduce further 
the system Eqs. (B.4) to five independent equations. 

We finally define the "la" uncertainty on the photon transfer function as £T 7 = 
[Rj — T^] 1 / 2 . However, the distribution of T 7 is not gaussian (and in general is also 
asymmetrical) so that T 7 ± 5T 7 should be interpreted just as a qualitative "band" with 
a not defined confidence level. 

For a constant value of a, the system of Eqs. (B.4) can be integrated analytically. 
For simplicity, we give here only the solution for the case a = 

_ . . 49 + 50e" 3 ^ 2 + Qe-'y 

&r(y) = Jog • (B.6) 

For y — > 00 we have thus the prediction T 7 = 2/3 ± l/3\/5. 
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